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Abstract. In this paper we consider construction of decoherence free subspaces for two- 
access random unitary channels. First, we concentrate on hermitian unitary noise model 
U for a bi-unitary channel and show that in this case a code exists if the space of Schmidt 
matrices of an eigensubspace of U exhibits certain properties of decomposability. Then, 
we show that our technique is also applicable in a generic case and consider its application 
two tri- (quatro-, . . . ) unitary channels. 



1. Introduction 

Quantum information transmission [TJ [21 [3l HI [5] through a non trivial 
quantum channel inevitably involves occurrence of errors. These errors, han- 
dled in an incompetent way, may completely shadow an intended quantum 
message. In this context, methods of combating such errors naturally emerge 
as one of the main topics of the theory of quantum channels. Luckily, several 
useful techniques have been developed to overcome destructive influence of 
coupling to the environment (see, e.g., [HE]). Among them, quantum error 
correction codes (QECC) are most widely recognized (see [8] and references 
therein). Methods of constructing QECC for quantum communication have 
been previously reported in the literature [9j [TUl [TT1 [12] . However, all these 
proposals concerned bipartite communication. No general approach has been 
developed to treat the case of the larger number of users of a quantum net- 
work. 

The present work is a step in an effort towards one such general ap- 
proach. We concentrate on random unitary two-access channels and show 
how to construct decoherence free subspaces (DFS) for such channels. DFS 
are the carriers of quantum information on which it is completely safe from 
the influence of an environment, that is a quantum message goes undisturbed 
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through a channel (in the sense of a solely unitary operation) . Since the recog- 
nition of their significance, DFS have drawn much attention (see, e.g., [13] 
and references therein) and the concept found its experimental realizations 
[HI [151 [16]. In our case, we have to protect two signals from spatially sep- 
arated independent senders, thus our DFS have to be factorizable according 
to this separation. It turns out that in the considered setting existence of 
DFS is exactly equivalent to the certain decomposability properties in the 
space of matrices of bounded rank. One of the motivations for the present 
paper was merging these concepts in a fruitful way. 

The paper is organized as follows. In Section 2. we provide some back- 
ground material. In Section 3. we state our main results. Next, we apply the 
results to the specific case of a qutrit-qutrit input. Then we discuss some 
generalizations. 

2. Background 

In this section we give some necessary background material. This includes 
quantum error correction conditions, spaces of matrices of bounded rank, and 
several miscellaneous mathematical facts. For the reader's convenience this 
part is quite voluminous. 

2.1. Quantum channels 

Quantum channel £ is a completely positive trace-preserving map. Every 
channel admits the so-called Kraus (or Choi-Kraus, or operator-sum) repre- 
sentation as follows C{g) = Y2i -AiQ-Ai with Yli^i^i = ^ [EMU]- A random 
unitary channel is the one which has the representation C{g) = ^2iPiUigUj , 
where Ui are unitary and YliPi = Pi — 0- When such a channel has two 
Kraus operators, i.e., C(g) = pUigUi* + (1— p]UigUy, *, it is called bi-unitary. 
This kind of channels are the main interest of the present paper. 

Channels can be classified upon the number of senders and receivers using 
them. We have the following types of channels according to such a classifi- 
cation [1(13 [201 EEl [22]: 

• bipartite — one sender and a single receiver, 

• multiple access — several senders and one receiver, 

• broadcast — one sender and several receivers, 

• fcm-user — k senders and m receivers (k,m > 1). 

In our reasonings we concentrate on two-access channels, that is multiple 
access channels with two senders. 
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Due to the possibility of a global rotation u\C{q)U\ on the output of 
a channel, both in bipartite and multiple access case one can consider a 
simplified bi-unitary channel in general reasonings^ 

£(q)=pq+(i-p)UqUK (1) 

For two-access channels it holds g = g\ (g) Qii where qi is an input of the i— th 
sender. 



2.2. Quantum error correction 

QECC is a subspace C of a larger Hilbert space %. Equivalently, a code 
is defined to be the projection Pq onto C C %. One says that C is correctable 
if all states g = Pc qPc can be recovered after action of a channel using some 
decoding operation V, that is VoC(g) = g. Such recovery operation exists if 
and only if PcA\AjPc = ctijPc for some hermitian matrix [L\%j = ay. These 
conditions are due to Knill and Laflamme (KL) [23J. 

When we have a larger number of senders we talk about local codes Cj 
prescribed for every sender. It is an immediate observation that KL condi- 
tions need only a little adjustment to serve for the case of MACs. Namely, 
we have (with the obvious notation): 

OBSERVATION 1. Local codes C{ arc correctable for a MAC with Kraus 
operators {A{\ with k inputs if and only if 

Pc x ® Pc 2 ® ■ ■ • ® P Ck A\A } P Cl ® P C2 ® . . . ® P Ck = aijP Cl ®Pc 2 ®...®Pc k (2) 
for some hermitian matrix [L]ij = ay. 

This is true since the set of product codes is a subset of the set of all codes. 
In further parts, we use the denotation R <g) R 1 or S (g> S' for a code for a two- 
access channel and talk about M ® N codes, where M, N denote dimensions 
of local codes. 

In case of many usages of the channel, A{ are tensor products of Kraus 
operators in KL conditions. In this paper, however, we concentrate on a 
single usage of a channel. For one use of a bi-unitary channel, Eq. ([1]), KL 
conditions reduce to a single condition PUP = \P which for MACs takes 
the form [H] 

R (8) R'UR ®R' = \R® R'. (3) 

In some cases there is no need to perform recovery operation as the output 
of the channel is already the same as the input. In these cases, code spaces are 



1 Naturally, in a concrete case one needs to remember that U = UvUz. 
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called decoherence free subspaces [13] . It is known that for bi-unitary channels 
DFS correspond to eigenspaces of degenerate eigenvalues of U [25]. Due to 
our interest in two-access channels we wish to find DFSs which factorize 
according to spatial separation of senders 0. 

It is worth explicitly noting that our problem is exactly equivalent to find- 
ing a product subspace in a given subspace. One may think of this problem 
as of the generalization of the well known (solved) problem of existence of a 
product state in a given subspace [26], [27]. 

2.3. Spaces of matrices of bounded rank 

A space of matrices of bounded (equal) rank is a space which contains 
only elements whose ranks are bounded by some prescribed number (or, 
besides the zero element, equal to it). The research on such spaces dates 
back to works by Flanders [28] and Westwick [29] . In quantum information 
theory the concept of spaces of matrices of bounded rank were recently used 
in Ref. [30]. If a space contains only elements of rank equal to k (besides 
zero element) we will be talking about A;-spaces. 

We will use 21, 23, £, . . . to denote space of matrices spanned respectively 
by Aj, Bi,C{, — We define r m (2t) = maxMgar(M). Further, we recall the 
concept of equivalence of spaces of matrices [28J and decomposability |29j 
We say that 21 is equivalent to 03 if there exist nonsingular matrices E,F 
such that 21 = {EBF, B £ 93}. If these matrices are explicitly specified 
we talk about (E,F) -equivalency. A subspace 21 of a x b matrices is called 
(t, s) -decomposable if it is equivalent to a subspace whose all elements have 
the form 



When t, s are not specified, we will be just talking about a general fact of 
(t + s) -decomposability. Below we collect several important facts concerning 
decomposability in the spaces of 3 x 3 matrices of bounded rank. It holds 
true that [231I3I1G3]: 

• three dimensional 2-subspaces are not 2-decomposable; moreover, up 
to an equivalence, there exists the unique such space — the space of 
skew-symmetric matrices H, 

2 It is obvious that if we need no correction at the output, each of the local subspaces 
was send noiselessly. In this context, a question arises whether it would make any sense to 
introduce the notion of local decoherence free subspaces. 

3 It is an interesting coincidence: this space corresponds to the antisymmetric space of 
two qutrits via the identification states-matrices; antisymmetric spaces play an important 
role in many applications of quantum information theory. 
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• if 25 is a 2-subspace then it necessarily holds that dim 33 < 3, 

• if a subspace 05 with r m (Q5) = 2 contains a rank one matrix then it is 
2-decomp osable , 

• a subspace 03 with r m (Q5) = 2 and 4 < dim's < 6 is 2-decomposable 
(follows from above). 

2.4. Space of states vs. space of Schmidt matrices 

One of our main tools will be the well known one-to-one identification 
of pure states with matrices: with fixed orthonormal bases \i) and \j) for 
Ui = C dl and H 2 = C d2 respectively, one defines the Schmidt matrix of a 
state l^) = c ijVj) to be C = J2ij c ijV)(j\- For two states \<f>) and \ip) with 
Schmidt matrices C and D it holds = trC^D. By direct calculation one 

finds that the transformation — > t/i <8> C^l^) corresponds to C — > UiCU^f. 
Further, we define the Schmidt rank r of a state \tp), denoted by r(|?/>)), to 
be the number of nonzero elements in its Schmidt decomposition and the 
maximal Schmidt rank r m of the subspace H to be r m (H) = max|^) GW r(\ip)). 
Obviously r(\ip}) = r(C). 

Let now % = span{|7j)} and = span{/ij}, where hi is a Schmidt matrix 
of 7j. Then r m (7i) = r m {$)) due to the isomorphism between the space of 
unnormalized states and matrices. It is thus natural to transfer the idea of 
decomposability of the space of matrices to the space of states. Therefore we 
propose to use the following 

PROPOSITION 2. A subspace % is called (i, j) -decomposable if S) is so. 



In our reasonings we consider channels with d dimensional inputs (and 
thus d 2 dimensional output). We take U to be hermitian, which implies that 
it must have the form 



where P and Q are both projectiono We further denote the eigensubspaces 
of U as V = PH and Q = QH and let dimP = p and dim Q = q. 

We have the following^ concerning the channel from Eq. (pQ) with U taken 
as in Eq. ([5]). 

THEOREM 3. A M ® N DFS exists if and only if at least one of the sub- 
spaces V and Q is (d — M, d — N) -decomposable. 

1 We naturally focus only on cases when U is not equal to 1. 

2 Note that a DFS for a bipartite channel with the considered noise model always exist. 



3. 



Main results 



U = P-Q 



(5) 
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PROOF. We put A = ±1 in Eq. (J3|) as these are the suspicious values (see 
Section [2X]) . We have using Eq. fl3J): 

R® R'{P -Q)R® R' = R® R' (A = 1) (6) 

and 

S®S'(P-Q)R®R' = -S®S' (A = -l), (7) 

where R, S and R' , S' are M and N dimensional projections respectively. It is 
evident that it is sufficient to conduct calculation for only one of the equations 
above as the sign of A exchanges only the roles of P and Q. Therefore we 
concentrate on Eq. ([6]), which due to P + Q = la can be rewritten in two 
equivalent forms as 

R ® R'QR ® R' = (8) 

and 

R® R'PR® R' = R® R'. (9) 

Let us concentrate on the first equation (second equation will be further in 
some cases more convenient to draw conclusions about existence of codes). 
Assume 



Q = Y f \<k){<k\, \<Pi) = Y, c kl\ kl )> k,i = i, — ,d, i = i,---,g (io) 



i=l kl 

and let [Ci]ki = c£j. We can represent the projections as rotated projections 
written in the standard basis 

Af-l N-l 
R®R' = u\®UlR®R!Ui®U 2 , R=Y,\9)(9\, & = J2\h){h\, (11) 

9=0 h=0 

where U\ and XJ% are unitary. Inserting Eqs (jlOD and ()lip into Eq. ([8]) and 
taking into account that a matrix is zero iff all its elements are so we arrive 
at 

(ui ® U 2 (j^ \<t> m ){<t>m\ \ u\ ® UU \kl) = 0, 
i,k = 0,l,...,M-l j,l = 0,1,..., AT- 1. (12) 
In particular, this must be true for ij = kl, which gives 

W m >| 2 = o, 

% = 0,1,...,M-1; i = 0,l,...,A-l (13) 
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and consequently for all values of m 



(ij\U 1 ®U 2 \<f> m )=0, 
i = 0,l,...,M-l; j = 0,l,...,N-l. 



(14) 



This condition implies that off-diagonal terms vanish as well and thus Eq. 
(TT4l) is fully equivalent to Eq. ([8]). Recalling the transformation rule for 
Schmidt matrices under local unitary rotations of a state (see Section 12.4. p 
condition above can be rewritten as 



Denote U = U\ and V = U-f- Eq. (I15p states that all Schmidt matrices d 
must brought with the same unitary U, V to a form with the M x N zero 
matrix in the upper left corner. Notice that instead of full unitary matrices 
U, V we can consider their reductions, i.e., isometries Ui som and Vi som , which 
are M x d and d x N respectively. Further, we need the following simple 
lemma 

LEMMA 4. Let A4 be complex d® d matrices and let V\, V2 be isometries. 
The condition V\AiV-2 = [0]mxAT holds for all values of an index 'i ' if and only 
if for all complex a = (a±, • • •) it holds that ViQ^ ajCj)V2 = [Q]AfxiV- 

This means that the space £ must be (U, V r )-equivalent to a space whose 
all elements have the zero M X N in the upper left corner. It suffices now 
to show that this equivalency means also (d — M, d — A r )-decomposability, 
i.e., to show that unitary matrices are as powerful as general full rank ma- 
trices in definition of decomposability (Eq. ©)• This what we are go- 
ing to demonstrate now. Let X, Y be non singular matrices. If XC{Y 
has a M x N zero block it means that X = ([X t ]^ x m, [X'] ( ix(d-M)) T an d 
^ = (IXldxJV) [Y'}dx(d-N)) where X and Y are respectively rank M and 
matrices such that XCiY = [0]mxN- It means that XCi has rank less than 
or equal to d — N and the (possibly nonorthogonal) rows of Y span the re- 
maining N dimensions allowed in a d dimensional space. Thus it is enough 
to take isometry whose rows span the same space as the ones of Y do to 
achieve zeroing of the resulting matrix. Arguing in the same way for left 
multiplication we arrive at conclusion that both X and Y can be replaced by 
isometries which can be completed to unitaries. This ends the proof of the 
theorem. ■ 

Let now D be the space spanned by Schmidt matrices of the projection 
P. We can immediately conclude the following simple necessary condition 
for the existence of a DFS 



i = 0,1 



{i\U 1 C m Uj\j)=Q,\/ m 
...,M-l;i = 0,l,...,A-l. 



(15) 
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COROLLARY 5. If there is a M®N DFC then either of the following holds 

r m (€) < 2d- (M + N), (16) 
r m (2>) < 2d- (M + N). (17) 

Alternatively 



r m (Q) < 2d -(M + N), 
r m (V) < 2d -(M + N). 



(18) 
(19) 



PROOF. Follows from Theorem [3]] and the property r(B) < t + u — r(A) — 
r(C) + r(ABC) valid for arbitrary matrices A sxt , B txu , C uxv . The latter 
can be proved by using twice the Sylvester's inequality |34|.B 

For further practical purposes, it is convenient to rewrite Theorem \3] as 
follows 

COROLLARY 6. A M (g> N DFS exists if and only if either of the following 
holds: 

i) there exist isometry [V®]mx<i such that 



( 



\ 



I V?C X \ 

v?c 2 
\ v?c q J 



\ 



<d-N, 



(20) 



ii) there exists isometry [Vf]Mxd such that 



( 



V 



VfD 2 
V VfD P J 



\ 



< d — N. 



(21) 



We observe that we can give immediately a rough bound on dimensions 
of an input for which a M ® N code surely exists for a fixed q. 

COROLLARY 7. Ifi)q = 2 and d > M + N or ii) q > 3 and d > min{Mg + 
N, Nq+M} then there always exists a M®N DFS for the considered channel. 

PROOF. The first part follows from the generalized Schur decomposition 
theorem |35j stating that for any two complex matrices A, B there exist 
unitary transformations U, V such that UAV and UBV are both upper or 
lower triangular at the same time. 

As to the second part. M rows of q Schmidt matrices of Q span at 
most Mq dimensional subspace. Since d > Mq + N there is still enough 
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space for N orthogonal vectors in the whole d dimensional space. Taking 
Uisom = (iMx(d-M) [0]mxm) and V isom to consist of the mentioned N 
vectors as matrix columns we produce M x N zero matrix. Analogous rea- 
soning applies to right multiplication by Visom and choosing Ui som to consist 
of proper columns. ■ 

Naturally, the same holds for p in place of q. 



We now move to the specific case of a qutrit inputs to the channel. Our 
main interest will be in the construction of 2 (8 2 codes. The cases 1 (8) 2 
(2(g) 1) are not interesting for us as such codes correspond to a non-zero rate 
only on the single wire of the network, on the other hand codes 2 (8 3 (3 (8 2) 
correspond to the situations when one of the parties can send with maximum 
rate, which clearly requires very specific type of noise (although not of the 
form U' g) 1 or 1 ® V as one can check) . 

Direct application of Theorem [3] relying on finding proper isometries 
Uisom and Visom giving zero blocks in the corner of all the matrices is a 
tedious task when 3 < q < 7. In these cases, one is thus recommended to 
consult Refs. [291 EH S3] f° r elegant techniques. Some relevant results have 
been given in Section 12.3.1 

In what follows we concentrate on the q = 2 case. In Ref. |31j . it has 
been shown that every two dimensional spaces of matrices with r m < k is 
^-decomposable. However, no explicit distinction between different kinds of 
decomposability has been given. We thus feel inclined to provide detailed 
discussion concerning (1, l)-decomposability in the relevant case of d = 3. 
The result is summarized in the following theorem. 

THEOREM 8. If q = 2 then a 2(8 2 code exists if and only if r m (C) < 2 and 
£ is not a (0,2) or (2,0) -decomposable 2-subspace. 

PROOF. (The proof is constructive) It is obvious from the upcoming Lemma 
I4T1 that we need only to solve Eq. ([6]), that is to take A = +1. Thus, no 
considerations concerning decomposability of 2) are required. Further, it is 
also evident from Corollary [3] that it must be that r m (<£) < 2 implying 
that r(Cj) < 2. We assume that we have already passed to the (W\,W2)- 
equivalent space where W\ and W2 come from the singular value decomposi- 
tion of C2. In such basis we assume these matrices to be 



4. Applications: qutrit-qutrit input case 




(22) 
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It is further assumed that r(Ci) < r{C 2 ). 

Suppose r m (<£) = 1. It is a simple observation that in this case, in the 
matrix C\ only the last row or column (not at the same time obviously) is 
non zero. Existence of a code is thus a trivial fact. 

Suppose now that £ contains an element of rank 2. W.l.o.g. we can 
assume that C 2 is such an element, that is both a and b are greater than zero. 
Condition r TO (£) < 2 is fulfilled if for all /3 it holds that det(d + f3C 2 ) = 0, 
which happens when0 



en = 0, det C\ = 0, oci 2 c 2 i + bc\ 3 c 3 \ = 0, 



(23) 



which we hereafter assume to hold. We now further consider specific cases. 
Suppose (ci2,ci 3 ) / (0,0) and (c 2 i,c 3 i) ^ (0,0). Take 




V, 



Q 



( 1 









QC21 

N 2 
bc 31 



\ 



V N 2 J 



(24) 



one can ver- 



12 



with normalization constants Ni = \J a 2 |ci2 1 2 + & 2 |ci 3 | 2 and N 2 
With this choice of isometries, using last condition of Eq. (|23p 
ify that it holds for the matrix elements that \V^C 1 V 2 Q ) n = [V^CiV^] 
[V^CiV^] 2 i = and V^C 2 V 2 Q = [0} 2x2 ; using last two conditions of Eq. 
(|23|) one arrives at [V®CiV®] 22 = 0. This proves existence of a code. 

If (ci2,ci 3 ) = (0,0) and (c2i,c 3 i) = (0,0) by previously mentioned gener- 
alized Schur decomposition and a simple swapping of rows we can transform 
both matrices simultaneously according to 




which outputs the proper code. 

Suppose (ci2,ci 3 ) = (0,0) and (c 21 ,c 31 ) / (0,0). Setting [Vf 
have 

/ Vl 2 C 2 l + Ul 3 C 3 l V 12C22 + Vl 3 C 32 V l 2 C 23 + V 1 3 C 33 \ 

C Q = V 22 C 2 i + V 23 C 3 i V 22 C 22 + f23C 3 2 V 22 C 23 + Ul 3 C 33 

v bvi 2 avi 3 

\ 6-U22 av 23 J 



(25) 



: Vij We 



(26) 



1 Conditions of this type are standard in analyses of spaces matrices of bounded rank. 
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which must be rank one (all rows, equivalently columns, must be proportional 
to each other). There are two (in principle non exclusive) possibilities for this 
to hold: (i) bv±2 = avis = and bv22 = 0^23 = or (ii) vyic<i\ + V13C31 = 
and t>22C2i + V23C31 = 0. The first alternative cannot be true in any case since 
a > 0, b > and this would entail the fact that vu = 173 = t>22 = ^23 = 
which is impossible because of the isometric character of V® . We thus have 
W12C21 + W13C31 = and V22C21 +W23C31 = 0. W.l.o.g. we can set v 2 2 = V23 = 
reducing the problem to finding conditions under which there exists such 7 
that the system of equations 

V12C21 + V13C31 = 0, 

jbvi2 = V 12 C 22 + V 13 C 32 , (27) 
jaVi 3 = V12C23 + -W13C33, 

has a nontrivial solution, where the last two equations are the requirement 
that the first row is proportional to the third one. This is possible if there 
exists such 7 that 

r(( 021 C22 " 7 ° C23 A)<1. (28) 

V V C 31 c 32 C33 - 76 J J 

Recalling that the first row of C\ is now equal to zero, it can be equivalently 
written as r{C\ — 7C2) < 1, which is simply the obligation for the subspace 
Q to contain a rank one element (there must exist such 7, not for all of 
them it must hold) or, in other words, not to be a 2-subspace. The case of 
(ci2,ci3) 7^ (0,0) and (021,031) = (0,0) can be treated obviously in a similar 
manner and we get that there must exist such 7 that 

Cl2 C13 

C22-7& c 23 (29) 
C32 C33 - 76 

with the same conclusion. This exhausts all possibilities. ■ 

With no effort we can extend our result to the case q = 7. Namely, we 
have 

THEOREM 9. If q = 7 then a 2 <g> 2 code exists if and only if r m (2>) < 2 
and D is not a (0,2) or (2,0) -decomposable 2-subspace. 

Notice that at no point in the proof have we made an assumption about 
orthogonality of the matrices. 

We conclude this section with some observations. 



FACT 10. Let P, Q, R be projectors with r(P) = p, r(Q) = q, r(R) = pq. 
The following holds P <g> QRP ®Q = P®Q if and only if R = P&Q. 
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This means that in some situations we can approach the problem of de- 
ciding existence of a code more directly for q = 4,5. Specifically, for q = 4 
we check whether Q is product 2 (8) 2. Positive answer resolves the matter on 
the spot. If the answer is negative, we need to check 2) for the decompos- 
ability . By analogy, if q = 5 we check whether P is 2 <S> 2 product, if it is we 
immediately have a code, if not — we test <£ for decomposability. Using this 
method we can easily show that A = — 1 and A = +1 cannot belong to h.2®2 
at the same time. We call this impossibility the uniqueness of DFS. We state 
this fact as follows: 

OBSERVATION 11. 2 <g> 2 DFS for the considered noise model in case of 
d = 3 is unique. 

Interestingly, the uniqueness is more powerful as the following holds true: 

OBSERVATION 12. If there is a 2 <g> 2 DFS for d = 3 there exists no code 
corresponding to other values of X. 

The proof of this fact is given in [37] . 

It should be noted that for the noise model introduced in the beginning, 
standard (non-product) DFS always exist, which follows from high degener- 
acy of spectrum. 

4.1. Examples 

Here we provide several illustrations to the results obtained above. We 
assume Eq. (jlOp to hold in what follows; in every case our concern is the 
existence of a 2 ® 2 code. We itemize our examples below. 

• q = 1 

|0 1 ) = -L(|OO) + |11) + |22)) (30) 

A code does not exist since the necessary condition (Corollary [3?]) is 
not fulfilled as r{C) = 3. 

• q = 2 

= -L(|ll> + |22» 

|&) = -^(|10> + |21» 

£ is a (0, 2)-decomposable 2-subspace and as such it is not (1,1)- 
decomposable thus a code does not exist (see Theorem |4"1) . 
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>2 



^=(|02> + |10» 
^(|01> + |20» 



£ is a 2-subspace but its (1, l)-decomposability can be easily seen. A 
code exists and is given by Pc = R R' = P\2 ® P12, where P\% = 
|1}(1| + |2)(2|. 



q = 3 

|0 1 ) = -L(|O1)-|1O}) 

|0 2 ) = -L(|12> - |21» 

|0 3 ) = ^(|02) - |20» 

Q is a projection onto the antisymmetric subspace. A DFS does not 
exist. This projection corresponds to the noise with U = SWAP. In 
fact, there is also no other code for two-access transmission through 
such a channel (see [37]). 

q = 4 

I'M = |oo) 
\<h) = |01) 
l^) = |io) 

|0 4 ) = |11) 

Q is product itself so a code definitely exists. Alternatively we could 
check £ and D which will result in (1, l)-decomposability of D. 

q = 4 

|0!) = -L(|00> + |H)) 

102) = |20) 
\<fo) = |21) 
104) = |22) 

It holds that r m (Q) = 3 so Q cannot be (1, l)-decomposable. On the 
other hand Q itself is not product so there is no code at all. It is 
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interesting that Q contains a rank three vector 1/V3(|00) + |11) + |22)) 
but its complement contains vectors of rank at most two [38] ■ The 
latter leads to 

OBSERVATION 13. In a 3® 3 Hilbert space, the complement V of the 
four dimensional subspace Q with r m {Q) = 3 may have r m (V) = 2. 

5. Generalization 

In this section we show that our reasoning can be applied to a more 
general type of noise. In fact, it turns out that a generic case in which a 
bipartite decoherence free subspace exists, which serves as a starting point of 
our discussions, can be treated by our technique. W.l.o.g. we focus on 2 ® 2 
codes in the d = 3 case. 

As we have already mentioned, a bipartite DFS exists only if the spectrum 
of U is sufficiently degenerate, in present case it must be at least four-fold 
degeneracy. This observation immediately leads to the conclusion that we 
can successfully consider the following noise model 

U = P + Y, ei8kp k (31) 

k^O 

with at least rank four projection Pq and arbitrary phases <5,. If at least three 
of them are distinct, it follows from the theory of higher rank numerical range 
that only A = +1 needs to be considered. If 5i = tt for i = 1, 2, 3, 4 then both 
A = — 1 and A = +1 have to be taken into account. Note that a phase in 
front of Pq is redundant as it corresponds to a global phase which cancels 
when U(-)U' is applied (this phase irrelevance also applies to the noise we 
considered so far). 

To support our claim concerning validity of the previous analyses, we 
provide a quick calculation . Let us assume that only A = +1 is possible. 
Since Yli -Pi = 1, we can replace U with U = 1 + J2k^o( elSk ~ ^)^k- Inserting 
this into the KL condition with A = 1 we obtain 

R ® R'(t + ( ei4 - VP^R ®R' = R®R' (32) 

k^k 

which further gives 

R ® R'{^2(e iSk - l)P k )R ® R' = 0. (33) 

Assume now that for each k we have P k = Ylm b^TH^fe 1 !' then for all a, (3 
it must hold that 

- l)<»^|C7i U 2 ^ mmut ® U}\afl = 0. (34) 

k m 
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Considering diagonal terms i.e., ij = a/3 leads to 

E( ei5fe - x ) E Ktfl^i ® ^ItDI 2 = 0. (35) 

k m 

Since Re(e l5fe — 1) < (phase equal to zero has been previously eliminated) 
we conclude that 

Vfc^Vij (^■|C/ 1 ®[/ 2 | 7 ^) = 0, (36) 

which, if we denote by C™ Schmidt matrices of j™, as previously is equivalent 
to 

Vfc^Vij (zlC/iC^li) = 0. (37) 

This means that the preceding analysis is valid for any U with at least four- 
fold degeneracy of the spectrum. 

We can take a step further in our generalizations and consider the triu- 
nitary channel 

A(gi <8> 02) = (1 - P - q)Qi <8> Q2 + pU (qi <g> Q 2 )U ] + qV{gi (g> g 2 )V ] (38) 
with 

U = P o + E ^"Q", V = P + ^2 e* k Qk 
k k 

with r(Po) > 4, aj / and arbitrary Qfc (orthogonal to Po)- This stems 
from the fact that in this case we need to solve system of equations stemming 
from KL conditions 

R R'UR ® R' = XiR ® R' (39) 
R R'VR <g> R' = \ 2 R ® Rl (40) 
R® R'U^VR® R' = \ 3 R® R'. (41) 

Taking A« = +1 we arrive at the same type of equations as previously. Due 
to the specific form of the unitaries, it is sufficient to consider only one of the 
equations and for the rest conclusion is drawn automatically. Increment in 
number of unitaries is naturally possible. 

6. Conclusions 

Summarizing, we have considered construction of decoherence free sub- 
spaces for biunitary two-access channels. Starting with a hermitian noise 
model we showed that the problem is directly related to the characterization 
of spaces of matrices of bounded rank, an area of linear algebra, which is 
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quite well established and understood. We further generalized the approach 
to arbitrary noise models and triunitary channels. To the best of our knowl- 
edge this the first time when analyses of spaces of matrices of bounded rank 
enter the field of quantum error correction. 

It is not clear to what extent the techniques we used can be applied to 
channels with arbitrary Kraus operators. Notice, however, that using our 
technique we can solve Eq. ([3]) for any operator with a spectrum with fixed 
real or imaginary part. 

We encourage the reader to consult Ref. [37J for methods of constructing 
higher entropy codes. 
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